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Abstract 

We discuss the bounds on polarized parton distributions which follow from 
their definition in terms of cross section asymmetries. We spell out how the 
bounds obtained in the naive parton model can be derived within perturbative 
QCD at leading order when all quark and gluon distributions are defined in terms 
of suitable physical processes. We specify a convenient physical definition for the 
polarized and unpolarized gluon distributions in terms of Higgs production from 
gluon fusion. We show that these bounds are modified by subleading corrections, 
and we determine them up to NLO. We examine the ensuing phenomenological 
implications, in particular in view of the determination of the polarized gluon 
distribution. 
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1. Introduction 

In the naive parton model, parton distributions are interpreted as probability densi- 
ties for partons to be found inside a given target [jT| . The partonic viewpoint immediately 
implies that if the individual parton distributions are bona fide probability densities, they 
should all be positive. Consequently, polarized distributions, which are differences of distri- 
butions of partons with polarization parallel or antiparallel to their parent hadron, should 
be bounded by the corresponding unpolarized distributions, which are sums over polar- 
izations. Because unpolarized distributions are known rather accurately, such bounds are 
potentially useful in the determination of the polarized distributions, and in particular of 
the polarized gluon distribution, which at present can only be indirectly extracted from 
scaling violations 0. It is well-known that, in perturbative QCD, the naive parton model 
is only recovered in the asymptotic limit, when all perturbative corrections can be ne- 
glected, while at finite energy the naive parton model results, in general, receive calculable 
perturbative corrections. It is thus interesting to ask to which extent the naive positivity 
bounds for polarized densities are modified by QCD corrections, especially in the case of 
the gluon distribution, which is of special phenomenological interest. 

Here we address this question by deriving bounds on the polarized quark and gluon 
distributions from positivity of cross sections. We show that the naive positivity bounds can 
be derived from leading order (LO) QCD, provided that the quark and gluon distributions 
are defined in terms of polarized and unpolarized physical processes. While it is clear 
that an adequate choice of defining process for quark densities is in terms of deep-inelastic 
structure functions, such as, e.g., F2 and gi, the selection of a suitable process for the gluon 
is subtler. In fact, in order to directly obtain the naive bound on \Ag\/g at LO, we need 
a process where only gluons of a given polarization contribute at LO. Thus, one is led to 
consider the production of a scalar particle (such as the Higgs) in gluon-proton collisions: 
only the gluon in the proton contributes at LO, and, moreover, if the external gluon and 
the target are polarized, only one polarization state for the gluon parton is selected. 

Next, we discuss how these bounds get modified by next-to-leading order (NLO) 
corrections, in general in a scheme-dependent way. In a generic scheme (such as the 
commonly used MS scheme) the NLO bounds may be more or less restrictive than the 
LO ones, depending on the parton distribution and the region of x which are considered. 
We then determine the full set of NLO bounds for quark and gluon partons explicitly. 
This requires the determination of the NLO corrections in the polarized case to the Higgs 
production process which defines the gluon at LO. After having established the improved 
NLO bounds, we use them to extract phenomenological information on polarized parton 
distributions in terms of the unpolarized ones, and specifically to constrain the polarized 
gluon distribution. 

2. Positivity of parton distributions at leading order 

In perturbative QCD, positivity bounds on parton distributions follow from their defi- 
nitions in terms of cross-section asymmetries, and from positivity of physical cross sections. 
Let us first consider the simple case of the structure function (71 (x, Q^) for inclusive deep- 
inelastic scattering of longitudinally polarized leptons on a longitudinally polarized target. 
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which leads to bounds on polarized quark and antiquark distributions. Neglecting power 
corrections, gi is related to the asymmetry Ai for deep-inelastic scattering of transversely 
polarized virtual photons on a longitudinally polarized nucleon through P] 



- <^3/2 _ 9lix,Q'^ 



cri/2 + Cr3/2 Fi{x,Q'^) 



(2.1) 



where the subscripts denote the total angular momentum of the photon-nucleon pair along 
the incoming lepton's direction. This immediately implies that gi is bounded by its unpo- 
larized counterpart Fi: 

|^7i(x,Q')| <Fi(x,g2). (2.2) 

The structure functions Fi and gi are related to parton distributions through coeffi- 
cient functions according to 



Uf 

i=l 



(2.3) 
(2.4) 



1=1 



where denotes the usual convolution with respect to x, and qi are quark and antiquark 
distributions of flavour i and with electric charge e^, {e'^) = Yl^=i , fi* is the gluon 

distribution, Ag^, Aqi and A^f the corresponding polarized distributions, and AC"^ and 
are, respectively, polarized and unpolarized coefficient functions (the index dis a reminder 
that these are deep-inelastic coefficient functions). The coefficient functions can be written 
as series in with x-dependent coefficients: 



oo ^ 



(2.5) 



fc=0 



Now, at LO 



(0) 



(x) = ACf ^°^(x) = 5(l-x) 



C^^(o)(x)=AC^' (°)(x) = 0, 
so that at LO eq. ( p.2| ) immediately implies 



(2.6) 



J2el [Aq,{x,Q^) + Aq,{x,Q^)] 



<^e2[g,(x,Q2)+^^(^,Q2)]_ (2.7) 



i=l 



The bound eq. (|2.7|) must be satisfied as a matter of principle for any choice of target, i.e. 
for any combination of quark plus antiquark distributions, and thus it must be satisfied 
by each quark flavour separately. Furthermore, the bound must also be satisfied by the 
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structure functions for charged-current scattering, and thus by each quark and antiquark 
distribution separately. Therefore, we conclude that at LO 

\/^q^{x,Q'')\<q,{x,Q'') (2.8) 

for all flavours z, for all x, and for all such that the LO approximation makes sense. 

The constraint eq. (|2.7|) immediately implies that at LO it is possible to interpret the 
polarized and unpolarized quark distributions according to a naive parton picture. Namely, 
if we define distributions of partons polarized parallel or antiparallel to their parent hadron 
by 

gf = ^fe±AQ,), ?f = ^fe±AQO, (2.9) 

then eq. (|2.8|) implies that qf and qf" are positive-(semi) definite and can thus be inter- 
preted as probability distributions. 

Next, one may ask whether a LO positivity bound can likewise be derived for the 
gluon distribution. In order to do so, we must select a pair of polarized and unpolarized 
processes which at LO define the gluon, i.e. processes whose LO coefficient functions are 
in an analogous way given by 

Cf'^°)(x) = ACf'^°)(x)=0 

(2 10) 

C^'(°)(x) = AC^'(°)(x) =5(l-x). 

Because the gluon only couples strongly, the leading nontrivial order must be at least 
0(as), unlike the quark which can be defined at 0{a^g). 

Processes which are typically used to determine the gluon distribution experimentally 
by singling out the photon-gluon fusion diagram, such as heavy quark production, are not 
useful to our present purpose, because even though they do have the property that at LO 
only the gluon coefficient functions are nonvanishing, they do not satisfy the requirement 
that the spin-fiip coefficient function vanishes, i.e. Cg = ACg. Therefore, with this 
definition positivity of the physical cross sections does not automatically imply positivity 
of the parton distribution. One could also think of using deep-inelastic scattering mediated 
by longitudinally polarized photons, i.e. in the unpolarized case the structure function F^, 
which receives a LO gluon contribution (of order ctg) Q (of order a^. However, there is 
no suitable polarized counterpart of F^. 

Nevertheless, the bound 

|A(7(x,g2)| <^(2:,q2)_ (2.11) 

is directly obtained at LO if we define the gluon distribution from inclusive Higgs produc- 
tion in gluon-proton scattering, i.e. from the process g+p ^ H + X. At leading nontrivial 
order [i.e. 0(a^)] this process proceeds through g + g —>■ H, induced by coupling of the 
gluons to a top quark loop (see fig. la). Even though this process is not realistic for an 
experimental measure of the gluon distribution, it is particularly convenient to formally 
define the gluon distribution and establish positivity bounds. The cross section for the LO 
partonic unpolarized process is given by 

a(,!,^H) = "-^<^'7< \A\H(s^ml), (2.12) 



288^2- 



n 
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Figure 1: Diagrams for the processes (a) g + g ^ H, (b) g + q ^ H + q, (c, d) 
g + g^H + g. 

where s is the center-of-mass energy of the gg coUision, and ^ is a dimensionless function 
of the Higgs mass and the top quark mass: A — A{m^/m1), such that hm^^oo |^('')|^ = 1- 
Furthermore, since the Higgs particle is scalar, only the polarized amplitude in which the 
two gluons have the same helicity is nonzero at LO in any frame where the two gluon 
momenta are in opposite directions. 

We can thus define the polarized and unpolarized gluon distributions by considering 
the gluon-proton scattering process. It is important to note that the external gluon is to 
be treated as a hadron, in the sense that it must be endowed with its own nonperturbative 
parton densities gg, Agg, Ug, Aug and so forth. The coUinear singularities associated with 
the initial gluon leg have then to be absorbed into a redefinition of the associated parton 
distributions. We are however free to choose the resulting distributions in a form which 
is convenient for our conceptual purposes. Our choice corresponds to that of an effective 
free gluon, i.e. gg{x) = Agg{x) = S{1 — x), Ug = Aug = . . . = 0. As a consequence, the 
gluon-proton cross section can be written as 



a[gp^ H + X]{x,ml) 



288\/27r 



A|2 




9{y,ml) 




(2.13) 



Aa[gp^ H + X]{x,ml) 



288A/27r 




where x = ^ and the singlet quark distribution is defined as 

n/ 

E(x, Q2) = J2 Q') + q^{x, , (2.14) 



and hkewise for the polarized distribution AE. 

The unpolarized and polarized cross sections a and Aa are respectively the sum and 
the difference of the cross sections with fixed parallel and antiparallel proton and gluon 
helicity. We have included a contribution from the scattering of the incoming gluon off 
quark and antiquark partons in the proton, which is present at higher orders even though 
it vanishes at LO where only the diagram of fig. la contributes. If we then expand the 
coefficient functions (i.e. rescaled partonic cross sections) C'* and AC'* in powers of as it 
follows immediately that the LO term in the expansion is given by eq. (|2.10| ). Positivity 
of the cross sections with fixed gluon and proton helicities then implies 

\Aa{x,ml)\<a{x,ml) (2.15) 

and thus the LO bound eq. ( |2.11| ). We conclude that the naive partonic positivity bound 
eqs. (^T8| ) , ( pTTTI) hold for quarks and gluons provided the scale is large enough that NL and 
higher order corrections are negligible. 

Since parton distributions themselves are scale dependent, one is led to ask how the 
positivity bounds will vary with scale at LO 0. The scale dependence of parton distri- 
butions is expressed in terms of the singlet distribution eq. (|2.14|) and nonsinglet quark 
distributions defined as 

n f 

?Ns(a:, Q2) = J2{w)-^) ^^*^^' + ^*(^' (2-16) 
1=1 

and likewise for the polarized distributions Aqns- The moments of the nonsinglet distri- 
butions then evolve multiplicatively: 

9Ns(iV, Q') = ?Ns(iV, Ql). (2.17) 



where for a generic distribution f{x) we have defined f{N) = dxx^~^ f{x). Moments 
of the singlet quark evolve by mixing with the gluon, according to 

[ci{N)E{N,Q^)+cl{N)g{N,Q^)] = 

a ^02^^-±W (2.18) 
Uglj) [ciiN)m.Ql) + cl{N)g{N^Ql)] , 

and likewise for polarized distributions. 

Now, the LO anomalous dimensions (which are the only relevant ones asymptotically) 
satisfy a+{N) > aNs{N) > a-{N) as it is trivial to verify from their explicit expression |jl|]. 
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It follows that the largest singlet anomalous dimension a+{N) will asymptotically drive 
the evolution of all the quark distributions for each flavour and of the gluon distribution. 
The same of course applies to polarized anomalous dimensions Aa{N). It is also easy to 
see that necessarily aj^{N) > |Aa"'"(A^)| for all A^: this is an immediate consequence of the 
fact that LO unpolarized and polarized anomalous dimensions are respectively sums and 
differences of the square moduli of invariant particle emission vertices (up to virtual 
contributions which are the same in the polarized and unpolarized case). Indeed, it is easy 
to check from their explicit expressions |jl|] that for all finite A^, aj^{N) > |Aa_|_(A^)|, while 
limjv^oo o^+(-^) = liniAf^oo ^ct+i^))- It follows that for all finite 

|A,.(Ar,Q^)| |A,(Ar,Q^)| 

Hence, the LO bounds eqs. (|2.(j|),( pTTTD will always be satisfied at large enough scale. 
Notice that, conversely, this implies that the LO bounds will always be violated at suffi- 
ciently low scale, or, otherwise stated, the (scale-independent) LO bounds become more 
and more restrictive as the scale is decreased, i.e. they become unacceptably restrictive if 
imposed at an excessively low scale. On the other hand, this raises the issue of computing 
subleading corrections to the bounds. We will do this in the next section. 



3. Positivity beyond leading order 

At NLO and beyond helicity-flip cross sections are no longer zero, so that the polarized 
and unpolarized coefficient functions C and AC are no longer equal to each other. In 
addition, the gluon DIS coefficient function Cg and the quark Higgs production coefficient 
function become nonzero. As a consequence, the LO positivity bounds, eqs. ( |2.8p ,( [2TTl 



receive NL corrections. In order to discuss these corrections, it is convenient to first rewrite 
the structure functions eqs. (|2.3|)-( p^ by separating the quark singlet and nonsinglet, 
rather than individual flavours, since the NLO bounds are affected by quark-gluon mixing, 
but this only involves the singlet quark distribution: 

Fi(x,Q2) = i^[cf,s®qNS + Cf(^E + 2nfC^®g], (3.1) 
gi{x, Q2) = i^[AC^5 ® ^QNS + ® AE + 2n/AC„^ ® Ag]. (3.2) 



The LO positivity bounds on individual quark and antiquark flavours eq. (|2.8| ) then 
imply LO positivity of the C-even combination 

|AQ,(a:,Q2)+AQ,(a:,g2)| < [q^{x , Q^) + q,{x , Q'')] , (3.3) 
which in turn implies LO positivity of the singlet quark distribution, eq. (|2.14|) 

|AE(x,Q2)| < i](x,Q2). (3.4) 

This is of course a necessary, but not sufficient condition for LO positivity of individual 
flavours, eq. (|3.3|), to hold. A necessary and sufficient condition for eq. ( p.3| ) is given 
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by constraining the polarized singlet by eq. ( |3.4|) , and the polarized nonsinglet through 
constraints which depend both on the value of AE and on the unpolarized (singlet and 
nonsinglet) distributions. The construction of such constraints is trivial: consider for in- 
stance the simple two-flavour case. The constraints eq. ( |3.3| ) then mean that a rectangular 
area is allowed in the (Au, Ad) plane. The corresponding equivalent constraints on AE 
and Aqisfs are found by giving coordinates in the plane with respect to the AJ^-Aq^s 
axes, which are rotated by 7r/4 with respect to the original axes. Likewise, for the separate 
positivity of quarks and antiquarks (eq. ( p.8|) ): necessary and sufficient conditions can 
be derived by imposing positivity of the C-even combination, eq. ( p.3|) , and a similarly 
derived condition on the C-odd combination Ag^ — Aqi. 

We are now ready to discuss the NLO modification of the quark and gluon bounds 
eqs. (p78|) , ( pTl] ) . First, note that the bounds on physical cross sections eqs. (p.2|)-( |2TT5D 
must be satisfied for all x, and therefore also by the corresponding Mellin moments: 

\gi{N,Q^)\<F^{N,Q^) (3.5) 
\Aa{N,Q^)\<aiN,Q^). (3.6) 

Since the bound eq. (|3.5|) must hold for any target, the same arguments which lead from 
eq. ( p.7|) to eq. ( |3.4|) now give the general bound on the singlet quark distribution 

\ACf{N, Q^)AE{N, Q2) + 2nfAC^{N, Q^)Ag{N, Q^)\ 

< [C^{N,Q')J:{N,Q')+2nfC^{N,Q^)g{N,Q'''' ^^'^^ 

which at NLO implies 

AS(iV, Q2) (i + ^ACt^ (^)(iV, Q2) j + 2nf^AC'g' ^'\n, Q')Ag{N, Q 



m, Q^) (l + tC-"' ^^\n, Q')) + 2nf^Ct (iV, Q^)g{N, Q^) 
Similarly, eq. ( |3.6| ) gives the NLO bound on the gluon distribution 

Ag{N,Q^) (l + fjACj'(^)(iV,Q2)) +^AC.'^'(^)(iV,g2)AE(iV,Q2; 



< 1 (3.^ 



< 1 (3.9) 



Equivalently, we may reabsorb the NL modification of the bounds in a change of 
factorization scheme. Indeed, we are always free to choose to work in a "parton" scheme in 
which the coefficient functions for the defining process, computed to any given perturbative 
order, coincide with their LO expressions 

Ct{N) = C;{N) = ACt{N)=ACi{N) = l 

(3 10) 

C^{N) = C^{N) = AC^{N) = AC^{N) = 0. 

Imposing that all higher order terms vanish C"^' ^"^^ = C^' = AC"^' = AC^' = 0, 
i > 1 determines the factorization scheme completely: given coefficient functions up to. 
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say, NLO in any scheme (such as the MS scheme) we can transform to the parton scheme 
by redefining the singlet quark and gluon distributions according to 



1 + 



a 



27r 



(jd, (1) 



(N) 
(AT) 



(N) 



h,{i) 



(3.11) 



and hkewise for polarized parton distributions. Otherwise stated, we are always free to 
choose a scheme in which the bounds reduce to their naive LO form. However, in a 
generic scheme parton distributions at NLO and beyond will not necessarily respect the LO 
positivity constraints eqs. (|2 
bounds eq. 



, (|2.11| ). Nevertheless, once we impose the correct positivity 
]9D on NLO parton distributions in any scheme of our choice, physical 
cross sections will be positive up to NNLO corrections, and if the parton distributions are 
transformed to any other scheme they will respect the positivity bound pertinent to that 
scheme. 

The coefficient functions C"^, AC^ and are all known up to NLO. In order to deter- 
mine the full set of NLO positivity bounds and study their phenomenological implications 
we must thus still determine the NLO polarized quark and gluon coefficient functions AC'*. 
We will do this in the next section. 



4. Determination of the polarized Higgs production coefficient functions to 
next-to-leading order 

At NLO the process g + p H + X can also proceed through gluon-quark fusion, 
g+q H+q (fig. lb), so the quark coefficient functions Cg , AC^ are nonzero; furthermore 
since there is an extra particle in the final state both helicity configurations can now 
contribute and the polarized and unpolarized coefficient functions are unequal. Similarly, 
the gluon-gluon fusion contribution can now proceed with emission of an extra gluon, 
g + g H + g (figs. lc,d), so that the gluon polarized and unpolarized coefficient functions 
are likewise unequal. The LO process finally receives virtual corrections which necessarily 
have the same momentum and spin dependence of the LO, and cancel IR divergences of the 
real emission contributions. All these contributions take a simple form when calculated in 
the limit ^ — > oo, in which the quark loop behaves as a pointlike effective interaction ||^. 
We will determine the polarized cross section at NLO within this approximation, which, 
regardless of its phenomenological relevance, is adequate for conceptual purposes. The 
corresponding unpolarized calculations have been performed in ref. [Of for the real emission 



diagrams while the virtual corrections have been determined in ref. [|T^ . 

The polarized parton cross sections which define the coefficient functions eq. ( |2.13| ) 
are determined by computing helicity differences of squared matrix elements, i.e., for the 
emission process g + f ^ H + f , with f = q, g we have 

AMf{s,t,u) = \My{s,t,u)\'^ - |M)^(s,t,w)p, (4.1) 

where s, t and u are the Mandelstam invariants of the process and the arrows refer to the 
helicities of the incoming g and / partons. This difference can in turn be determined by 
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contracting the tensor amplitudes with polarized projectors (more precisely, differences of 
projectors): for the real gluon emission process g + g ^ H + g 

AM,{s,t,u) = -M^^^M;-P,P^^^P^^p, (4.2) 

where the gluon projector is given by 

P^u — '''^tivpcr ^ _ ^ (4-3) 

in terms of the momentum k carried by the gluon line on which the projector is applied 
and an arbitrary lightlike vector n such that k ■ n ^ (for example, the momentum of the 
other gluon line) . For the quark emission process g + q H + q we have 

AM,(s, t,u) = Tr [Mj^P'i ^ MJ- P^^, (4.4) 

where the trace runs over fermion indices, k and k' are the incoming and outgoing quark 
momenta, and the quark projector on the incoming line is given by 

P, = -75- (4.5) 
We can now compute AMf in dimensional regularization in d dimensions; the calcula- 



tion closely parallels that of the unpolarized contribution P-p!0| . The polarized calculation 



requires a recipe to handle the 75 matrix and the antisymmetric tensor: we will continue 



them to d dimensions according to the HVBM [|TT| prescription, since this will give results 
in the same scheme in which the NLO anomalous dimensions of all parton distributions 
have been determined In the gluon case the tensor amplitudes are given in ref. 0, 
while in the quark case we can readily calculate the amplitude from the diagram of 
fig. lb. We get 

AM,(M.«)^«2G^.5^!<±4^ (4.6) 

fi4 — 1/^ 

AM,(., t, u) = ^iGf^'-^, (4.7) 

where Gp is the (i-dimensional continuation of the strength of the effective ggH interac- 
tion [ITO. The form of this continuation is immaterial because it appears as an overall 



factor in all the separately divergent contributions which add up to the final finite NLO 
result. 

The kinematics for the inclusive process can be entirely specified in terms of m|, 

2 

s = and the center-of-mass scattering angle of the emitted parton cos^ = 1 — 2y, in 
terms of which the remaining Mandelstam invariants are given by 

t = -{s~ ml){l -y) = -ml{l - y) 

1 - X ^^-^^ 
u = -{s- m^)y = -ruf^y . 

X 
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The physical cross-section is obtained by integrating the real emission amplitude over phase 
space, adding vitual corrections and collinear and ultraviolet counterterms, and supplying 
the appropriate flux factor: 

Aa((7 + / ^ + /) = ^ y"#2 AM/ + AaY' + Aaf + Aaf . (4.9) 

The two-body phase space in (i = 2(2 — e) dimensions is given by 



1 1 f Atx\^ f ^ 



l-2e 



The collinear counterterms remove the singularities which correspond to collinear 
radiation of the flnal-state parton by the incoming one, and are thus given by 

Aaf'= ['^Aao{mly)-T{l + e){4nrP,f(-) (4.11) 

Jx y e \yJ 

where Pgf is the LO QCD splitting function 0, the factor of r(l + e){4TvY accompany- 
ing the e pole characterizes the MS subtraction scheme, and Aao{mf^,x = rn\/ s) is the 
polarized Born cross section in d dimensions, which turns out to coincide with the unpolar- 
ized four-dimensional cross-section eq. (|2.12|) up to the replacement Gp ^ G^.lll Virtual 



corrections, as well as the ultraviolet counterterm related to charge renormalization, only 
contribute to the gluon process, since the quark process only starts at NLO, and are pro- 



portional to the Born cross section in d dimensions. They were computed in ref. from 
which they can be simply obtained by replacing the d dimensional unpolarized Born cross 
section with its polarized counterpart, and are thus given by 



Aa^* = r(l + e){AnYAao{ml x)[^-^-^-2n^ 



(4.12) 

ll-2/3n/ 1 
Itt 2e" 



Aa7 = Aao(m|, x) ^f^-r(l + e)(47r) 



Performing the phase-space integral eq. ( |4.10| ) of the square amplitudes eq. ([4.7|)-( ^r6| 



and adding the collinear counterterms ( [4.11 ) and, in the gluon case, the virtual corrections 
eq. (|4.12| ) we get the physical cross sections for the NLO processes g + g H + g and 
g + q H + q. Double poles in e cancel between real and virtual contributions to the 
gluon process. Poles in e proportional to 5(1 — x) are generated in the real contribution to 
the gluon process because of the presence of a term proportional to (1 — x)~^^^'^^'^ which 
has an IR singularity in the e — * limit: 

' - '-S{l-x) + -^-2e(^^] +0{e% (4.13) 



;i-x)i+2^ 2e ' ' (l-x)+ \ 1-x 



^ Notice that the fact that the dependence on e of Actq only comes through G% is nontrivial, and 
peculiar of the polarized calculation where it follows from the adoption of the HVMB prescription 

for the antysimmetric tensor. 
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These poles cancel against those in the virtual correction eq. ( 4.12 ). The collinear simple 
poles are cancelled by the splitting function counterterms in both the quark and gluon 
emission diagrams. 

The final finite result for the physical cross sections gives us immediately the coefficient 
functions AC'^ eq. ( p.l3| ) by dividing out the physical Born cross section eq. ( |2.12|) : 



11- 



:i 



X 



+ (11 + 27r2) 5(1 - 



X) 



+ 2 



ln(l 



X 



X 



ARgg{x) — \nX APgg{x) 



+ 21n(^)AP,,(x) 



X 



,(1 



X] 



X 



+ APgg{x) 



In 



X 



+ In 



X 



where APfg are LO QCD polarized splitting functions, 



ARgg{x) 



:i 



X. 



+ 



APggix) - (^S{1 - X) 



ml 



(4.14) 
(4.15) 



(4.16) 



/3o = 11 — |n/ is the leading coefficient of the QCD beta function, and fx^ is a factorization 
scale. We can now derive the full set of positivity bounds for the quark and gluon dis- 
tribution. The result eq. ( [4.14|) - (|4.15|) is of course also in principle interesting for its own 
sake as a determination of the cross section for Higgs production in polarized gluon-proton 
scattering.! 



5. Phenomenology of next-to-leading order positivity bounds 



We can now collect all results for the eight coefficient functions which determine the 
NLO positivity bounds for quarks and gluons, eqs. (|3.8| ) -( ^^ , in moment space. The Fi 
coefficient functions MS can be determined from the known F2 [13| and Fl coefficient 
functions using Fi = {F2 — Fl)/(2x). The gi coefficient functions are available fl^ in the 



MS scheme with HVBM prescription; from these it is straightforward to determine them in 



the AB version of the MS scheme [jT5| which allows a comparison to recent determinations 
of parton distributions f^. H The full Higgs production unpolarized coefficient functions 



^ Note that eqs. ( [4.14D -( p:.15P give the coefficient of ^ in the perturbative expansion of the 
coefficient function. The coefficient functions eqs. (4.14)-( ^.15| ) are thus related to the 0{^^^ 



a.. 



contributions to the cross section for the process g + f^H + X hy Aa{g + f H + X) 



^ The AB scheme is a variant of the MS-HVBM scheme in polarized DIS which adopts a 
different subtraction in the definition of the polarized gluon coefficient function. The scheme is 
constructed in terms of the MS scheme in such a way that only the Zqg entry of the scheme change 
matrix is nonzero. This means that only the ACg coefficient function is affected by the scheme 

change, while all other coefficient fucntions are the same as in the MS scheme. 
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Figure 2: The eight NLO coefficient functions, evaluated at Q^=l GeV'^, as functions 
of the moment variable N. 



are given in ref. , while the calculation of their polarized counterparts was presented in 
the previous section. For completeness, we list in Appendix A the full set of NLO polarized 
and unpolarized coefficient functions in moment space in the MS-AB scheme. 

Combined bounds on A'E{N, Q"^) and Ag{N, Q"^) can now be derived using eqs. (|3.8| )- 
( p.9|) for various scales and values of the moment variable. We will assume everywhere 
the equality of the factorization and renormalization scales; specifically in eq. ( 4. 14 ), ( 4. 15 ) 



/i = m^. The bounds are only significant for > 1, because the first moments of 
unpolarized distributions diverge. Furthermore, the scale has to be chosen large enough 
that the NL truncation of perturbation theory is a reasonable approximation. On the 
other hand, the bounds will only be relevant at not too large values of the scale, since 
at asymptotically large scales all bounds are trivially satisfied because of the asymptotic 
vanishing, eq. (|2.19D , of the ratios ^ and Furthermore, at asymptotically large Q^, 
NLO corrections are negligible and all bounds reduce to the naive LO ones discussed in 
Sect. 2. 

In fig. 2 the eight NLO coefficient functions are plotted for = 1 GeV^ as functions 
of the moment variable. It is interesting to note that the NLO corrections to the Higgs 
production process in the gluon channel (i.e. to AC^ and Cg) are very large (as large as 
100% at the Higgs scale, and much larger at the low scale considered here). As pointed 
out in ref. [|10|, this is mostly due to the very large size of virtual corrections. These 
however, are proportional to the Born process and thus only contribute to 



corrections, 

the gluon channel spin-nonfiip cross section, and indeed. 
Since the NLO corrections to the positivity bounds are only due to the nonvanishing of 



' g 
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the spin-flip process, it follows that the NLO corrections to the bounds are actually very 
small and can be safely treated within perturbation theory despite the large size of the 
NLO corrections to the coefficient functions themselves. 

The pattern of deviations from the LO behavior displayed by the NLO corrections to 
the coefficient functions of fig. 2 is rather complicated: all the Higgs production coefficient 
functions satisfy \/\C^\ < for small N but \AC'^\ > when is large; the DIS 
coefficient functions for all have |ACf | < but |AC^| > (notice that both ACg 
and Cg are negative). The fact that some NLO corrections to the polarized coefficient 
functions are larger than their unpolarized counterparts means that some of the NLO 
contributions to the associate cross sections are negative. 

The bounds eqs. (p78|)-( |3T9D are useful to constrain more poorly known parton distri- 
butions by the better known ones. Specifically, the bounds are significant for high values of 
the moment variable because this singles out the large x region, where the unpolarized 
distribution is reasonably constrained by the data whereas very little is known about the 
polarized distribution. 

In fig. 3 we present the LO and NLO bounds eqs. (|3.8| )-( |0|) for several relevant values 
of the moment variable N . The bounds are given in the form of an allowed region in the 
AE-A^f plane. The contours depend on the unpolarized distributions, which are taken 



from the CTEQ4LQ set . The maximal allowed values of the unpolarized distributions 
are thus affected by the uncertainity on the unpolarized distributions. The size of the 
uncertainty of parton distributions cannot be evaluated precisely at present. However, in 
the central x region (0.01 ^ a; ^ 0.8) the uncertainty can be considered to be essentially 
negligible on the unpolarized quark distribution, while a recent analysis suggests that 



the uncertainty on the unpolarized gluon distribution is of order 10% unless x is very small 
(x ~ 10"'^) or large {x ~ 0.2) the gluon being vanishingly small at large x. This should 
give an idea on the expected uncertainty on the bounds of fig. 3 due to the errors on 
parton distributions, while the difference between LO and NLO bounds gives an idea of 
the uncertainty related to higher order corrections. 

The bounds are presented at the scale = 1 GeV^ which is often chosen as "initial" 
scale in parametrization of parton distributions, and which can be taken as a lower limit 
for the applicability of perturbation theory. The NLO corrections are seen to modify and 
distort the allowed area in the AS- A^r plane. The pattern of modification of the LO bounds 
follows the pattern of NLO corrections to the coefficient functions displayed in fig. 2. In 
particular, it is clear that if |AC(^)(A^)| < C(^)(A^) (i.e. the NLO contributions to the 
various cross sections are all positive) then the NLO bound is less restrictive than the LO 
bound. For instance the quark bound becomes 

[C'^{N, Q')E{N, Q^) + 2nf {C^{N, Q^)g{N, Q^) + AC^^(7V, Q^) A^(iV, Q^))] 

<AE(Ar,g2)< (5.1) 

[CfjN, Q^)S(Ar, Q^) + 2nf {C^jN, Q^)g{N, Q^) - AC^jN, Q^)Ag{N, Q^))] 

ACi{N,Q^) 
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Figure 3: The LO (dashed hnes) and NLO (soHd Unes) positivity bounds on AT,{N) and 
Ag{N) for = 1 GeV^ and N = 1.5, 2, 5. The values of AS and Ag corresponding 
to the NLO fits to gi data of ref. p| are also shown. 
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This is always less restrictive than 



|AS(iV,Q^)|< ^^^^^^^^^^ S(iV,Q^) (5.2) 

because C^(iV, Q'^)g{N, Q'^)±AC^{N, Q'^)Ag{N, Q^) is always positive if \AC^\ < since 
to NLO accuracy we can use eq. for the gluon contribution. Eq. (|5.2|) is manifestly 

less restrictive than the LO bound eq. (U) if \ACf\ < . If instead \ACf\ > the NLO 
bound will be more restrictive if the gluon contribution can be neglected (so that eq. ( |5.2| ) 
holds) but could still be more restrictive than the LO one if the gluon contribution is large. 
Similar arguments apply to the gluon distribution. 

Inspection of fig. 3 shows that at = 1 GeV^ the NLO corrections to the LO bounds 
turn out to be reasonably small, but non-negligible for small A^-moments, where they tend 
to make the bound less restrictive than the LO bound if both AS and Ag are positive (as 
experiments appear to indicate) . The corrections are instead essentially negligible for large 
values of A^. However, at sufficiently low scale all the NLO corrections will start growing 
and dramatically modify the LO bounds, but use of perturbation theory at such scales is 
not advisable. 

Figure 3 also shows the values of AS and Ag obtained in a recent determination 0] of 
polarized parton distributions. Different fits correspond to different choices of functional 
forms and theoretical assumptions on the parton distributions, and give a feeling for the 
spread of the results. It is apparent that one fits (fit B) display large percentage violations 
of the positivity bound on Ag{N) for large values of the moment variable A^. This reflects 
the fact that large values of correspond to the large x region where Ag{N) is extremely 
small and the fitted results are affected by very large uncertainties. Indeed, the main 
purpose of ref. was instead to provide a reliable determination of the first moments of 
AS and Ag, while the large x i.e. large A^ behavior of the parton distributions could not 
be determined accurately. This however shows explicitly how positivity bounds can be 
used to improve a determination of parton distributions, and specifically Ag by providing 
independent information. Indeed, we could modify fit B in such a way that the bounds are 
respected, without significantly affecting the quality of the fit to data and the determination 
of the first moment. However, imposing that the bound is satisfied we would get a fit with 
reduced uncertainity on A^r in the large x region. 



6. Summary and outlook 

In the naive parton model the polarized parton densities in a proton with helicity +1/2, 
and p_, are interpreted as probability densities and, as a consequence, are positive 
semidefinite. Thus, given Ap = — p_ and p = p^ + p-, they satisfy \Ap\ < p, for any 
p = Qi, Qi, g. The QCD-improved parton model reproduces the naive parton results at LO, 
i.e. when all higher order perturbative corrections (as well as all power corrections) are 
neglected. In this limit, if the naive positivity bounds for parton densities were violated, 
one could exhibit hard processes with unacceptable negative cross-sections. For quarks 
obvious examples of such processes are given by the polarized asymmetries Ai = gi/Fi 
eq. (|2.1| ). For gluons, the selection of the process which leads to the bound \Ag\ < g is less 
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direct. One needs a process whose rate on a polarized proton with hehcity +1/2 starts 
with Qj^ and does not contain either or q± or q±. The process g-\-p H + X, where H 
is a neutral, colourless massive scalar (for example, a Higgs boson) satisfies the required 
properties. 

As is well known, most naive parton model properties and sum rules are in general 
violated if higher order QCD corrections are included. For example, the charge sum rules 
u{l,Q'^) — 'u(l,(5^) = 2 and similar ones [where u{N,Q'^) is the N-th moment of the u- 
quark density u{x, Q"^)] are in general violated, unless the unpolarized quark densities are 



defined by the structure function F2,[|T^ i.e. by extending to all orders, by definition, the 
LO expression of F2 in terms of quark densities. Similarly also the momentum sum rule 
or the LO positivity conditions are in general violated at NLO and higher accuracy, with 
corrections which depend on the definitions of the relevant quantities beyond LO. 

In this article we have computed the NLO corrections to the positivity relations in 
a set of specified definitions. In particular we adopted the MS scheme and a simple set 
of physical processes for the definition of parton densities, including the above mentioned 
Higgs production process for the gluon definition. Of particular interest in practice are the 
positivity bounds in the polarized singlet sector AE, A^r. This is because the determination 
of the polarized gluon density attracts considerable attention at present. When fitting the 
data on scaling violations often the gluon comes close to the LO positivity bound at large x, 
so that the precise form of the bound at NLO can be relevant. While at LO the positivity 
bounds in the singlet sector are given separately for AE(A^) and Ag{N) (|AE(A^)| < E(A^), 
|A5f(iV)| < g{N)), at NLO due to quark-gluon mixing the allowed domain is modified into 
a more complicated region of the AE — Ag plane. 

Our detailed quantitative results show that the modified region in the AE — Ag 
plane can be either more or less restrictive than the LO positivity domain, depending on 
the values of N, E(A^) and g{N). Interestingly, we also find that the corrections to the 
positivity limiting contours are in general much smaller than the size of the corrections 
on the individual coefficient functions. In particular this is true for the corrections to 
the Higgs process. It is well known that the NLO corrections to Higgs production by 
gluon-gluon fusion are quite large |1^. But the bulk of the corrective terms is from virtual 
diagrams which, being proportional to the Born term, do not modify the LO bounds. 

As a consequence, we find that, even at relatively small values of Q^, such as 
~ 1 GeV^, the modified domains are reasonably close to the corresponding LO re- 
gion. Of course if is further lowered down to a clearly non-perturbative region, the 
NLO corrections eventually explode and the NLO positivity contours become very dif- 
ferent than the LO contours. It is important to note that if the LO positivity bounds 
are imposed at very small they in general produce exaggerately constraining limits on 
the parton densities. This is because the evolution goes in the direction of making the 
positivity bounds trivial ai ^ 00. Hence, applying the LO bounds at very small 
leads to an evolved domain at larger that is considerably smaller than the LO bound 
at the larger value of Q^. It follows that if parton distributions are parametrized at a low 
scale, positivity bounds can only be safely imposed at a higher scale, where perturbation 
theory is reliable. A possible strategy is to use positivity bounds to discard fitted parton 
distributions which do not satisfy them, thereby reducing the uncertainty on the results 
of the fit. 
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Imposing positivity bounds consistently at NLO in a given scheme will then guaran- 
tee positivity of physical cross sections for the defining processes. It will also guarantee 
positivity of any other process whose LO coefficient functions do not have the partonic 
form eqs. (|2.6|),( prT0| ), because the LO bounds for such processes will be less restrictive. It 
will not, however, guarantee positivity of other processes for which eqs. (p.6| ), (p.lO|) hold. 



because then positivity would depend on the relative sign and size of NLO corrections. 

In conclusion, positivity bounds on parton distributions derived from specific defin- 
ing processes are in general useful as necessary conditions which provide complementary 
information on polarized parton distributions, especially at large x. 
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Appendix A. Coefficient functions in moment space 

We list here the MeUin moments C{N) = x^~^C{x) dx of the NLO contributions 
to the coefficient functions in the MS scheme for generic values of the ratio of factorization 

2 

and renormalization scales k = The full coefficient functions are defined in eq. ( |2.5|) ; 

Mj? 

the LO contributions are given in eqs. (p^ ) , ( pTTO|) . 
Deep-inelastic scattering NLO coefficient functions: 



S,iN-l)[S.iN + l) + l)+±-^^ 



^2(iV-l) 



Ci^^^\N) = -TR 



N{N + 1){N + 2) N + 2 



+ Pqg In/C 



ACs''^^\n) = Ci 



Si{N-l) \^- + S,{N + r 
N -1 



+ APqq{N) Ink 



ACf (i)(iV) = -Th-^-—^ [S,{N - 1) + 1] + AP,,(iV) \nk. 



9 
2 



(A.l) 

(A.2) 

(A.3) 
(A.4) 



The AB-MS scheme |T5[ (used in sect. 5) is obtained adding to ACg'^^^(A^) eq. a 
contribution equal to — while all other coefficient functions are left unchanged. 

Higgs production NLO coefficient functions: 



C^'(^\N)=Ca 



c 



h,il) 



11 + 47r^ 



22 



(N) = C 



3 N{N^ -1){N + 2)_ 

+ 2Qgg{N)+2Pgg{N) lli k 

- - 3 

'^77(iv^^ + ^^'^^^^ + ^^^^^^ 



11 + 47r^ 



+ 



22 



3 ■ N{N^ -1){N + 2)_ 

+ 2AQgg{N)+2APgg{N) lu/C 

^ +AQg,{N) + APgg{N)lnk. 



We have defined 



^ 1 



(A.5) 
(A.6) 

(A.7) 
(A.8) 

(A.9) 



k=i 
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Qgg{N) = 2Ca [Sl{N + 2) - 2Sl{N + 1) 

+ZSl{N) - 2Sl{N - 1) + Sl{N - 2)] 
/^Qgg{N) = 2Ca [2Sf{N + 1) - SSfiN) + 2SiiN - 1)] 



Qgg{N) = Cf 



45i (iV - 1) ^ 4:Si (AT) 2Si {N + 1) 



N -1 
2 



N 



N + 1 



+ 



1 



(AT - 1)2 Ar2 
45i(Ar) ^ 25i(iV + 1 



(iV+l)2 
2 



TV 



+ 



1 



N + 1 7V2 (AT + 1)2 



(A.IO) 
(A.ll) 



(A.12) 
(A.13) 



Furthermore, Pij{N) and APiJN) are the usual LO polarized and unpolarized anoma- 
lous dimensions, i.e. the Mellin transforms of the LO QCD splitting functions, whose 
expressions we list here for completeness: 



P — 



19 



P — 

^gq ~ 



AR 



Q9 



AP, 



91 



Cf 
Tr 

Cf 



3 _ ^ _ 1 
2~N~ N + 1 

+ N + 2 
N{N +1){N + 2) 

N'^ + N + 2 



-2Si{N-l) 



P = 

^99 



APg,= 



2Ca 
Cf 
Tr 

C F 



N{N^ - 1) 
1 



+ 



1 



N{N-1) (iV + l)(+2) 
3 1 1 
2~N~ N + 1 
N-1 



-Si{N) 



2 



2Si{N-l) 



(A.14) 



APgg = 2Ca 



N{N+1) 

N + 2 
N{N + 1) 
2 



N(N + 1) 



Si{N) 



2 ' 
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